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We have examined the validity of the Gor'kov expansion in the strength of the order parameter 
of type II superconductors near the upper critical field. Although the degeneracy of the electron 
levels in a magnetic field gives non-perturbative terms in the solution to the Bogoliubov-de Gennes 
equations we find, contrary to recent claims, that these non-perturbative terms cancel in the expres- 
sion for the thermodynamic potential and that the traditional Gor'kov theory is correct sufficiently 
close to H c 2 at finite temperature. We have derived conditions for the validity of the Gor'kov theory 
which essentially state that the change in the quasiparticle energies as compared to the normal state 
energies cannot be too large compared to the temperature. 



Recently the view has been advocated, that the Gor'kov expansion describing type II superconductors close to the 
upper critical field H C 2 may be invalidUl. Even within the mean-field approximation there is a possibility of non- 
perturbative effects arising from the degeneracy of the Landau levels for electrons moving in a magnetic field. This 
degeneracy means that even a small perturbation (e.g. superconducting order) can change the quasiparticle levels 
significantly as compared to the normal state. This effect has been proposed as a mechanism for the breakdown 
of the standard perturbation theory describing type II superconductors even close to H c i- This non-perturbative 
effect should give rise to effects such as tails of residual superconductivity above the usual H c2 , the possibility of the 
superconducting transition being first order, and unusual behavior of the heat capacity, magnetisation etc. close to 
the phase boundaryEJ. In this work we have examined this possibility. We show that for T = there are indeed terms 
not contained in the Gor'kov expansion for the difference between the ground state energy of the mixed state and 
the normal state. This is in agreement with the results obtained by BahcallEJ and is of no surprise since the Gor'kov 
expansion is essentially a high temperature series. However, for T^Owe show that the non-perturbative terms in the 
difference Qs — f2 jv between the thermodynamic potential in the mixed state and the normal state vanish and indeed 
that the Gor'kov expansion is a convergent series for the superconducting order A(r) not too large. We thus prove 
incorrect the claim of Bahcall that there is a non-perturbative third order term in the expression for Q$ ~ We 
have derived some criteria for the convergence radius of the Gor'kov expansion in the order parameter. A comparison 
between the results of the Gor'kov expansion and a numerical solution of the corressponding Bogoliubov-de Gennes 
(BdG) equationso confirms our conclusions. In this .paper we work in two dimensions. We do not consider any 
fluctuation effects relevant to high T c superconductorsB. 



Within the mean-field approximation the mixed state of a type II superconductor is described by the solutions to 
the BdG-equations. In a constant magnetic field the order parameter forms an Abrikosov vortex lattice and it is 
convenient to use a set of single particle states <^jv,k characterized by the Landau level index N and a wavevector k 
in the Brillouin zone of the vortex lattice. In this basis the BdG-equations split up into a 2N x 2N secular matrix 
equation for each k, where N is the number of Landau levels participating in the pairing. In this basis the BdG 
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where u v Nk is the coefficient of (f>Nk for the Bogoliubov function w k (r) and v^ k is the coefficient of </>^_ k for the 
function v k (r), and £jy = {N + l/2)hcu c — /i. fi is the chemical potential. The off-diagonal elements F-^mm are: 

Funm = J drA(r)0(r)^ k 0(r)^_ k (2) 

and the order parameter is determined self-consistently as: 

A(r)=gJ2<(r)vl(ry(l-2f2) (3) 

where g is the coupling strength and / k = (1 + exp(_E^/ fcgT)) is the fermi function. We neglect any finite Zeeman 
splitting for simplicity. Due to the translational symmetry the order parameter A(r) is completely characterised by 
a finite set of parameters A^u. For notational simplicity we work in the lowest Landau level approximation (LLL) 
(i.e Aj^o = 0) in which the center-of-mass motion of the Cooper-pairs has the kinetic energy Hlo c /2 where u c is 
the cyclotron frequency. None of the conclusions in this paper are altered when this restriction is relaxed. When 
the chemical potential /x is at a Landau level (i.e nj = fi/huj c — 1/2 = integer ) we have, in the normal state, exact 
degeneracy between an electron state in the Landau level n/ + m and a hole state in the Landau level n/ — m. 
Likewise when the chemical potential is exactly in between two Landau levels such that rif = n + 1/2 there is 
degeneracy between an electron in a level nf + m + 1/2 and a hole in a level n/ — m — 1/2. When this is the case 
we expect the possible non-perturbative effects of a finite order-parameter to be strongest. We will show that the 
convergence radius for the Gor'kov equations is indeed smallest when nf is an integer. To examine the validity of the 
Gor'kov expansion it is convenient to use the following expressions for the difference — fijv in the thermodynamic 
potential between the mixed state and the normal state: 

il s -n N = - ( dr|A(r)| 2 -2fc i 3T^ln(cosh(/3£; A r k /2)) + 2fc B TL>^ln(cosh(/3^/2)) (4) 

^ Nk N 

Here D = is the number of k-vectors in the Brillouin zone, V is the volume, B is the magnetic field, and 

(3 = 1/kgT. 



III. ZERO TEMPERATURE 



To illustrate the origin of the non-perturbative effect for T = it is sufficient to examine the case when only one 
Landau level participates in thepairing and rif is an integer. In this case the positive energy solution to equation (jl]) 
is E nf k — \Fknfnr\- Equation (|4|) reduces to: 

E gS -E gN = - ( dr|A(r)| 2 - ^£„ /k (5) 

k 

Since \Fk njni \ oc Ao oc |A(r)| we see that we obtain a linear term in |A(r)| in equation (Q). This is a non-perturbative 
term since the Gor'kov expansion only contains even powers of the order parameter. This T = result is unaltered 
when we have many Landau levels participating in the pairing and it agrees with the result obtained by Bahcalltl. It 
is a trivial consequence of the fact that we have to take the T — > limit fcgT ln(2 cosh(/3i?jvk2)) — > -Ejvk/2 before we 
perturbatively expand the result in the size of the order parameter. 



IV. FINITE TEMPERATURE 



A. Quantum limit 



For finite temperature the situation is different. It is now possible to expand ln(2 cosh(/3i5/v k /2)l.in powers of the 
order parameter and then check if we obtain any non-perturbative terms, as proposed by Bahcallj. For notational 
simplicity we wi ll again do the calculation in the quantum limit when only one Landau level participates in the pairing. 
In section IV B we will treat the slight modifications in our result when more than one Landau level are within the 
pairing width. The quasiparticle energy is now i?„ k = \/? 2 + | F^ nn \ 2 . We need to expand ln(cosh(/3i? nk /2)) in 
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|-Fknn| 2 - Writing [3E n \^/2 = \/e 2 + z 2 where e = /3£„/2 and z = (3\Fk nn \/2 we are lead to consider the analytic 
properties of the function ln(cosh(Ve 2 + z 2 ))- The poles and branch cuts in the complex plane zgC determine the 
convergence radius r$ for a power series in z. A simple analysis gives: ro = \J t 2 + it 2 /4. The requirement for the 
convergence of a perturbation series for ln(2 cosh(/3i? jyk/2)) is then 

\F knn \ < Va+^ 2 (k B T) 2 (6) 

This requirement is most restrictive when the Landau level is at the chemical potential (£„ = 0). Wc then have 

£k < k B Tn (7) 

Furthermore, we see that there will appear only even powers of |-Fknn| m the series. This is true for general fj, (i.e also 
when £„ = 0). So we have ruled out any non-perturbative cubic term in the expression for f2<j — fijv thereby disproving 
earlier predictions based on a numerical analysisa. Doing the expansion and comparing with a standard expression for 
the thermodynamic potential based on Gor'kov's equationsQ we find (not surprisingly) that it reproduces the Gor'kov 
series term by term. The convergence of the Gor'kov expansion is determined be equation (^). It is now clear that 
the Gor'kov expansion is a high temperature series. So the break down of the theory for T = is of no surprise. For 
finite T we expect the Gor'kov series first become unreliable when there is a Landau level at the chemical potential. 
This is because the requirement in equation (^) is most restrictive when = and because the superconductivity 
and thereby the change in the quasiparticle energies (obtained by a self-consistent solution of equation ([!]) is enhanced 
when there is a Landau level at the chemical potential]. 



B. Several Landau levels 



The above conclusions are essentially unaltered when there is more than one Landau level participating in the 
pairing. We calculate the quasiparticle energies from equation (Q) perturbatively in F-^ nm using degenerate and 
non-degenrate perturbation theory. Then we expand ln(2 cosh(/3i?/vk/2)) in powers of the order parameter. The 
convergence radius for the series is again smallest when the chemical potential is at a Landau level. The only 
complication is that we obtain both even and odd powers of F^ nm in the expression for the quasiparticle energies. 
But the odd terms cancel in the expression for Vis — SIn due to the fact that there are two quasiparticle levels when 
=/= for which the odd powers in the expression for the energy have opposite signs. There is only one positive 
energy solution for the case £„, = though. However, the odd terms from this solution vanish in the expression for 
Sis — VI n due to the fact that d 2l+1 ln(cosh(a;))| x= o = where I is an integer. A long tedious calculation shows that 
we recover the standard terms in the Gor'kov series. A sufficient condition for the convergence of the Gor'kov series 
is 

E nk - £ n < min [2k B Tn, k B T^(3 2 £ 2 + (8) 

which has to hold for each quasiparticle level within the pairing region. So we expect that the Gor'kov theory breaks 
down when significant portions of the quasiparticle bands lie outside the regions defined in equation (|§|) . It should be 
noted that one cannot ignore the contribution from higher quasiparticle levels to Qs — Ojv (£n 0)- This is easily 
seen from equation (|j) since 

ln(cosh(/3(£ + 8E)/2)) - ln(cosh(/3£/2)) > Mposh[p{5E)/2)) i,8E>Q (9) 

So any treatment which focuses only on the quasiparticle level at the chemical potential will ignore important contri- 
butions to Vl s . Equation (^) can be transformed into the requirement: 

< |A(r)| 2 >=^f dr|A(r)| 2 < 2^]^(k B Tf (10) 

Based on an extensive numerical analysis, Norman have suggested a similar conditiqa. Since the Ginzburg- 

Landau equations are derived from the microscopic BCS theory using the Gor'kov expansion^ it would be of interest 
to restate the above criterium in terms of Ginzburg-Landau parameters. Doing this we obtain: 

Hr2 ~ H <^ 1/2 (7C(3)) 1/2 /3^(.-^)(|) 2 (11) 



H c (0) 
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Here £(x) is Riemann's Zeta function, k is the Ginzburg Landau parameter, Pa is the Abrikosov parameter, and 7 is 
Euler's constant. This restriction is always fulfilled for type II superconductors within the normal range of validity of 
the Ginzburg-Landau equations (i.e \T — T c \/T c <C 1). 

As an example of the breakdown of the perturbation series we have plotted the orderparameter Ao as a function 

of rif = 1/2. In figure 1 we have plotted both a numerical exact solution of the BdG-equations and the fourth 

order perturbative result using a method developed earlier!. We have chosen the parameters such that Ud/oJc = 5, 
g/fiioj 2 = 8.2 and krjT '/Jvjj c = 0.3 when rif = 12. As can be seen the perturbation theory agrees fairly well with the 
exact solution. The perturbative result tends to differ the most from the exact solution when the chemical potential is 
at a Landau level (n/ integer). This is in agreement with the above remarks. To illustrate the temperature dependence 
of the convergence radius of the Gor'kov series we have in figure 2 again plotted Ao as a function of rif for a very low 
temperature. As can be seen the perturbation series breaks down much earlier Ao ~ 1500 for this low temperature in 
agreement with equation (||) . In figure 3 we have plotted the lowest quasiparticle level along a high symmetry direction 
in k-space when rif — 12 for the parameters used in figure 2. The horizontal line gives the boundaries for — £ 
calculated from equation (||). There are large parts of the quasiparticle band in k-space lying outside the region of 
convergence of the perturbation expansion thereby explaining the observed discrepancies between perturbation theory 
and the exact numerical result. 



V. CONCLUSION 



In conclusion we have examined the recently debated validity of the Gor'kov expansion. The conclusion is that 
although the degeneracy of the normal state levels gives large effects on the quasiparticle wavefunctions even for 
weak superconducting order these effects cancel in the expression for the thermodynamic potential and the Gor'kov 
expansion is correct for finite temperature. We have therefore ruled out the possibilty of a non-perturbative third 
order term in the expression for the thermodynamic potential. The range of validity of the Gor'kov expansion is 
given by equation (|^) which shows that it is essentially a high temperature expansion. This requirement is always 
fulfilled within the range of validity of the Ginzburg-Landau equations leading to no inconsistencies. Furthermore we 
have the usual requirement \E n k — £ n | <C Tiu> c for perturbation theory to work. We expect the Gor'kov expansion 
the break down first when the chemical potential is at a Landau level. Our results are confirmed by a comparison 
between the results of an exact numerical solution to the BdG-equations and a Gor'kov series to fourth order in the 
order parameter. 
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Figure Captions 



Fig. 1: The order parameter Ao vs nf — fi/fkJ c — 1/2 calculated numerically (solid line) and perturbatively to fourth 
order in Ao (dotted line) for lj,i/lj c = 5, h J [2 = 8.2 and kBT/hto c = 0.3 when nj = 12. 

Fig. 2: The order parameter Ao vs nf — (i/huj c — 1/2 calculated numerically (solid line) and perturbatively to fourth 
order in Ao (dotted line) for LUd/io c = 5, h J p = 7.0 and ksT/fiuJc = 0.05 when n/ = 12. 

Fig. 3: The lowest quasiparticle band in units of 7kj c along a high symmetry line in k-space for n/ = 12 . The dashed 
line marks the boundary defined by equation (g). 
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